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We present a multi-band model for superconductivity at the metallic interface between insulating 
oxides LaAlC >3 and SrTiOs (001). Using a self-consistent Bogoliubov-de Gennes theory, formulated 
with the realistic bands at the interface, we investigate the spin-singlet and spin-triplet pairings in 
intra-band and inter-band channels. We find that the Rashba and atomic spin-orbit interactions 
at the interface induce singlet pairing in the inter-band channel and triplet pairing in both the 
intra-band and inter-band channels when the pairing amplitude in the singlet intra-band channel 
is finite. The gate-voltage variation of superconductivity is resolved in different pairing channels, 
compared with experimental results and found to match quite well. Interestingly, an enhancement 
of the superconducting transition temperature by external in-plane magnetic field is found revealing 
the existence of a hidden superconducting state above the observed one. As the interface is known 
to possess high level of inhomogeneity, we explore the role of non-magnetic disorder incorporating 
thermal phase fluctuations by using a Monte-Carlo method. We show that even after the transition 
to the non-superconducting phase, driven by temperature or magnetic field, the interface possesses 
localized Cooper pairs whose signature was observed in previous experiments. 

PACS numbers: 74.78.-w, 74.20.Rp, 74.40.-n, 61.43.Bn 


I. INTRODUCTION 

The discovery of superconductivity ( T c ~ 200 mK) 
at the interface 1,2 between perovskite band insulators 
LaA10 3 (LAO) and SrTiOa (STO) triggered a plethora 
of investigations 3-12 in the last few years due to its ex¬ 
otic nature arising primarily from the presence of com¬ 
peting ferromagnetism (Tcurie — 200 K) 13-15 , spin-orbit 
interaction (SOI) 16,17 and disorder 18-20 . Besides, the 
quasi-two dimensional electron gas (q2DEG) at the inter¬ 
face 21 exhibits intriguing novel properties such as metal- 
insulator transition 22-26 and ferroelectricity 2 1 . On top of 
that, the ability to control these properties by external 
electric field 28,29 added an extra-dimension to the nano¬ 
electronics industry 30,31 . 

The q2DEG is formed by an electronic transfer mech¬ 
anism in which half an electronic charge per unit cell is 
transferred to the interface to avoid a polar discontinu¬ 
ity 32,33 . The electrons are confined in a few Ti02 layers 
located within a region of about 10 nm thickness at the 
interface and occupy the t 2 g orbitals of Ti ions 34,35 . Den¬ 
sity functional theory (DFT) reveals that the d xy band is 
situated below the d yz , d zx bands by ~ 0.4 eV due to the 
confinement at the interface 36,37 . The spin-degeneracy of 
the bands is lifted by the inversion symmetry-breaking 
Rashba SOI and an atomic SOI 38,39 . Magnetotransport 
measurements infer the presence of two types of carriers 
with different mobilities and the high-mobility carriers 
have been predicted to be responsible for superconduc¬ 
tivity 40-42 . Micliaeli et al. 3 suggested that the system 
hosts the antagonistic ferromagnetic and superconduct¬ 
ing orders by favouring a disordered stabilized helical 
FFLO state induced by strong Rashba SOI. The com¬ 
plex nature of the coexisting phases 13,14,43 has naturally 
led to the predictions of unconventional superconductiv¬ 
ity 44 and different magnetic ground states. The micro¬ 


scopic understanding of the origin of superconductivity 
remained obscure until the recent convincing evidence of 
electron-phonon coupling, obtained using tunneling spec¬ 
troscopy 45 . It is, therefore, apparent that phonons play 
the dominant role in electron pairing in other STO-based 
superconductors such as doped STO 46 and X/STO (X 
= LaTiOa 47 , GdTiC> 3 48 , FeSe 49 ) interfaces. The super¬ 
conductivity at the LAO/STO interface is unique in the 
following aspects: (i) it appears at very low career con¬ 
centrations (~ 10 -13 cm -2 ) 28 , (ii) the transition temper¬ 
ature ( T c ) shows BKT-like behaviour 28 while the pairing- 
gap or the superfluid density follows BCS prediction: 
2A o/iksTg) ~ 3.4, where A 0 is the pairing gap at T = 0, 
ks is the Boltzmann constant and T g is the so called ’gap¬ 
closing temperature’ 50 , (iii) it coexists with inhomoge- 
nous ferromagnetic puddles of large moments (~ 0.4 fj.s 
per interface unit cell) 13 . The coexistence of the compet¬ 
ing orders, albeit in phase segregated regions 4,5,43 , gives 
rise to fascinating phenomena such as the enhancement 
of superconductivity by magnetic field’ 1 and the mag¬ 
netic field assisted transient superconductivity 52 leading 
to possible ’hidden order’. However, despite intensive 
previous studies, complete theoretical understanding of 
the nature of the multi-band superconductivity is lack¬ 
ing and necessitates a careful and thorough theoretical 
analysis. 

In the following, we use a three-orbital model for su¬ 
perconductivity to develop an understanding of the na¬ 
ture of superconductivity at the interface in the pres¬ 
ence of magnetic moments and try to shed light on the 
questions raised above. We study the spin-singlet and 
spin-triplet electron pairing in intra-band and inter-band 
channels. Using a self-consistent Bogoliubov-de Gennes 
(BdG) method, formulated with the realistic bands at 
the interface, we explore the mean-field phase-diagrams, 
the role of spin-orbit interactions and external magnetic 
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field on the electron pairing. It is found that the pairing 
in the singlet inter-band channel and triplet intra-band 
and inter-band channels are induced by the Rashba and 
atomic SOI when the pairing amplitude in the singlet 
intra-band channel is finite. Taking cue from the ex¬ 
perimental data, we incorporate the gate-voltage in our 
analysis, study the gate-voltage variation of the pairing 
amplitudes and plot the phase-diagram to compare with 
the experimental results. We study the behaviour of the 
pairing amplitudes in presence of an external in-plane 
magnetic field and find an enhancement of the supercon¬ 
ducting transition temperature when the magnetic field 
is applied along certain directions in the interface plane. 
The magnetic field enhancement of superconductivity has 
been observed experimentally 51 and arises because of the 
interplay between superconductivity and the competing 
ferromagnetism. It suggests a hidden superconducting 
phase above the superconducting transition temperature. 
Since the interface superconductivity is highly inhomoge- 
nous in nature and appears at very low career concentra¬ 
tion, the thermal phase-fluctuation becomes significant. 
We study the phase-transition from superconductor to 
a non-superconducting state, driven by temperature or 
perpendicular magnetic field, taking into account the 
thermal phase-fluctuation using a Monte-Carlo method 
and observe that there are localized Cooper-pairs in the 
non-superconducting phase. The presence of these local¬ 
ized Cooper-pairs has been confirmed in previous exper¬ 
iment 53 . 

The remainder of this paper is organized as follows. In 
Sec. II, we discuss about the origin of ferromagnetism and 
superconductivity and introduce the multi-orbital effec¬ 
tive Hamiltonian for the interface q2DEG. In Sec. Ill, we 
formulate the self-consistent equations for the supercon¬ 
ducting order parameters within the multi-orbital BdG 
framework. In Sec. IV, we present and discuss our re¬ 
sults obtained within the self-consistent BdG method. In 
Sec. V, we present of Monte-Carlo analysis of the thermal 
phase-fluctuation in superconductivity. The conclusions 
are briefly summarized in Sec. VI. 


II. EFFECTIVE THREE-ORBITAL MODEL OF 
SUPERCONDUCTIVITY 

In the following, we elaborate the microscopic mech¬ 
anisms of ferromagnetism, superconductivity and other 
ingredients of the interface q2DEG and establish the ef¬ 
fective Hamiltonian for the interface electrons. 

The electrons coming from the top LaA 103 layer to 
neutralize the polarization discontinuity at the interface, 
predominantly occupy the d xy orbitals of Ti ions in the 
terminating Ti 02 layer and establish quarter-filled d xy 
states. Because of large onsite Hubbard, and nearest- 
neighbour Coulomb repulsive interactions at the inter¬ 
face, all the electrons get localized at the interface sites 
and form a charge-ordered insulating ground state with 
a weak anti-ferromagnetic super-exchange coupling me¬ 


diated via the Oxygens 3,54 . Additional electrons, sup¬ 
plied by the application of the back-gate voltage or the 
Oxygen vacancies near the interface, will find the top 
Ti 02 layer as energetically unfavourable and go to the 
next Ti02 layer to occupy the t- 2 g orbitals of Ti ions. 
The electrons in the t 2 g orbitals in the next T 1 O 2 layer 
participate in conduction and exhibit superconductivity. 
Spectroscopic experiment 42 and DFT studies 38,39 show 
that, due to confinement at the interface, the d xy band is 
lower in energy at the T-point by ~ 0.4 eV than the quasi- 
one dimensional, relatively heavier d yz ,d zx bands. When 
the Fermi-level is tuned, the system encounters a Lifshitz 
transition at an electron concentration n c ~ 1.68 x 10 13 
cm -2 where the low-mobility d yz ,d zx electrons start get¬ 
ting occupied 42 as depicted in FIG. 1. An interesting fea- 
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FIG. 1. (Color online) (a) The band structure of the three 
tig orbitals in presence of the Rashba and atomic SOI (spec¬ 
trum of IKo + 34 aso + Jtnso)- There is a mixing of the d xy 
orbital (blue lines) with the d yz orbital (green lines) and the 
d zx orbital (red lines) due to the Rashba SOI. The dashed 
horizontal line, at fi ~ —1.04 eV, denotes a Lifshitz transition 
point at which the two upper bands start getting occupied, 
(b) the total density of states as a function of energy. The 
Lifshitz transition is reflected by the sharp jump in the den¬ 
sity of states near the transition point. (c)-(d)-(e) depicts the 
change in the Fermi-surface topology at energy below, near 
and above the Lifshitz transition. 


ture of the band-structure is that because of the atomic 
SOI, the orbital characters of the d xy band and the heav¬ 
ier d yz ,d zx bands get interchanged near the Lifshitz point 
and the splitting due to Rashba SOI is significant only 
near the degeneracy points. The spin-degeneracy in all 
the bands is lifted by the spin-orbit interactions, the split¬ 
ting being largest near the band-mixing points. 

The wave-functions of the itinerant electrons in the 
t 2 g orbitals in the Ti 02 layer below the interface are 
extended to the terminating Ti02 layer. Therefore, 
these electrons interact via a ferromagnetic exchange 
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with the localized moments leading to the in-plane fer¬ 
romagnetic order. The exchange interaction can be de¬ 
scribed by a Hamiltonian 'Kfm = Ei Ea I drJ a S(R i ) • 
s a (r)S(r — R,), where S(R;) is the spin operator at the 
local moment sites R*, s a (r) is the spin-density oper¬ 
ator in the itinerant orbital a (d xy , d yzi d zx ), J a is the 
strength of the exchange interaction. Since the itin¬ 
erant d yz ,d zx orbitals are orthogonal to the localized 
d xy orbital, J ( i r Jd yz ,d zj; 3 - This is in agreement 
with the spectroscopic studies which indicate d xy na¬ 
ture of the in-plane ferromagnetism 55 . When treated 
at the mean-field level, "Kfm will essentially be reduced 

tO K FM = J2k,a,a,<r'( h XOL® X )v,<?' c \aa Ckaa '' where h xa. 
are the Zeeman splitting amplitudes, corresponding to 
different orbitals of index a, along the in-plane direc¬ 
tion (taken to be along x axis) with h xa = h x i for 
d xy orbital and h xa = h x 2 for d yz , d zx orbitals (h x 1 > 
h x 2 ). The attractive interaction, mediated by electron- 
phonon coupling, for the three itinerant orbitals can be 

expressed as R sc = -g Y,k,k\ a ,f) c Lt c -W c - fc '^ Cfc '“t> 
where g is the strength of the pairwise electron-electron 
interaction. The spin-singlet and spin-triplet pair¬ 
ing amplitudes in the intra-band and inter-band chan¬ 
nels can be defined as A® g = —g (ck a ^C-k@i) and 
A Q/ 3 a<T = ~9 (ckaaC-kp*) respectively. Neglecting 

fluctuations beyond mean-field, the pairing term be¬ 
comes R SC = Efc,a,/ 3 ( A «/ 3 C Lt C -fe/ 3 t + A a/ 3 tt C Lt ct -fc/ 3 t + 
A a/3u c Li ct -fe/3t + h.c.). We have Aj, m = A^ m 
and the orbitals are indexed according to (a, b , c) = 

xy 5 dyz , d zx ). 


the following Hamiltonian: 

^RSO = 7 51 ( C laa C kbo C kc* ) 

k,cr 

( 0 2 z sin k x 2 isin/cy \ / Ckacr \ 

2i sin k x 0 0 Ckba ^ 

2 i sin k y 0 0 / \ Ck ^ ) 

where 7 = 20 meV is the strength of the Rashba SOI. It 
is interesting to note that the Rashba SOI, described by 
R rsOi is very different from what is usually observed in 
the 2DEG at semiconducting hetero-interfaces. 

The Oxygen vacancies, which are developed at the 
interface during the deposition process, are considered 
as indispensable parts of the interface q2DEG and 
have very significant role in ferromagnetism, supercon¬ 
ductivity and their coexistence 4,5 . We model these 
non-magnetic impurities as the local random shifts in 
the chemical potential and express by the Hamiltonian 
Mdis = E^ Eaa Vldc l d aa c ida<y, where i d denotes the 
defect sites of total number N d , V ld is a random poten¬ 
tial which varies within a range [—IT, IT] and IT is the 
strength of the disorder. The percentage defect concen¬ 
tration is given by rid = Nd/N 2 x 100. 

The total effective Hamiltonian for the interface elec¬ 
trons is, therefore, given by 

Re// = Rq + RaSO + RflSO + RfM + R SC + Rdis (3) 


Another significant feature of the interface q2DEG is 
the presence of the atomic and Rashba SOI which re¬ 
organize the spin and orbital degrees of freedom of the 
t 2 g electrons. The atomic SOI, described by the Hamil¬ 
tonian R ^so = A so l ■ s, appears because of the crystal 
field splitting of the atomic orbitals. In the t 2 g orbital ba¬ 
sis (Cfcaf, c kb f, Cfcct, Cfcaj., c k bi, c kc i), CKaso can be written 
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where A so = 19.3 meV is the strength of the atomic SOI. 
On the other hand, the Rashba SOI, which describes the 
broken inversion symmetry at the interface, is given by 


where R 0 = Efc, Q ,o-( e fca “ vWkaa 0 *™ describes the 
band dispersion of the electrons in the three t 2 g orbitals 
with e ka = — 2 fi(cos k x + cos k y ) — t 2 — R 3 cos k x cos k y , 
e k b = —ii(l + 2 cos k y ) — 2 t 2 COsk x — 2t^cosk v , e kc = 
—ti(l + 2 cos k x ) — 2t 2 cos k y — 2t^ cos k x , /i is the chemi¬ 
cal potential and t\ = 0.277 eV, O = 0.031 eV, 1 3 = 0.076 
eV are the tight-binding parameters 38 . 

Although multi-band superconductivity in this inter¬ 
face q2DEG was proposed earlier 7,56 , the explicit nature 
of the pairing symmetry and the intra-band or inter-band 
superconductivity were not explored. Our model uses the 
realistic band structures, obtained from DFT studies 38 , 
and treats the electron pairing in intra-band and inter¬ 
band channels within the multi-band BdG theory, to be 
described below. 


III. MULTI-BAND BDG THEORY 


The self-consistent BdG theory is perhaps the best 
available numerical technique to study the interplay of 
superconductivity with real-space inhomogeneity or com¬ 
peting orders such as ferromagnetism or charge density 
wave within mean-field approximation in any experimen¬ 
tally realizable geometry. To begin with, the Hamilto- 
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nian ^K e f f in Eq. (3) is written in real lattice as 

K B dG = ~ C \** C 3P°' + h - c •) 

ij, a, (3, a,a' 

— ^2 (g — V d Sii d )cj arT Ci ar j 


E ( h * 


2,Q!,Cr,<7 


+ E ( A a/3(^)4at C W + 

2 , 0 !,/3 

+ E ( A L/3 fft r(n)cLa C ^a + h - C ') 

<ij> : C£,(3 ,a 


(4) 


where is the tight-binding hopping amplitudes 

which contains 1K 0 , and ^Krso , A®^ = 

-g (Ciai-Cipi) and = -g (c iQO -c i/3(7 ) are the local 

singlet and triplet pairing gaps, 5u d is the Kronecker’s 
delta function. 

The Hamiltonian ‘KsdG in Eq- (4) is diagonal¬ 
ized by the unitary Bogoliubov transformation c.i aa = 
U nL„ln + <aa7n f wllich J^ds the multi-band 
BdG equations (ex', being a dummy index, is omitted 
hereafter.) 
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(5) 


where - [{g - V id 8 iid )5 a(T > - 

(/i xa cr x ) (Ta ./](5,j(5 Q: /3. For a square lattice of size N x N, the 
BdG Hamiltonian matrix has dimension 12iV 2 x 12iV 2 for 
three orbitals. Using the above Bogoliubov transforma¬ 
tion, the local pairing gaps can be obtained as 


A a /3 = -f *fli - ul„,v % *a+] tanh 
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\2k B Tj 

( 6 ) 


where f(E n ) = 1/(1 + exp^jj/fegT)) is the fermi func¬ 
tion. In what follows, the BdG equations (5) are solved 
numerically on a square lattice with periodic boundary 
conditions to find out the eigenvalues E n and the local 
quasi-particle amplitudes u l nacr , u l rlarj and the new pair¬ 
ing gaps are calculated using Eq. (6). This process is 
repeated until self-consistency is reached at every lat¬ 
tice sites. Finally, the average values are obtained via 

A aP = B/ N ) £* A^. 

It is important to mention here that static mean- 
field theory neglects fluctuations and, therefore, over¬ 
estimates the fermionic field amplitudes. The temper¬ 
ature or the magnetic field, being treated in the analysis, 
are therefore shouldbe seen as just parameters and their 
qualitative features, not quantitative estimates, are rele¬ 
vant. 


IV. RESULTS 

Having formulated the self-consistent BdG theory, we 
now present the numerical results obtained on a 25 x 25 
square lattice. First we analyze, in the homogeneous sit¬ 
uation (rid = 0), the effects of Rashba and atomic SOI 
on superconductivity, the gap-structures and phase di¬ 
agrams of the singlet and triplet pairings and then the 
effects of in-plane magnetic field on superconductivity. 

Spin orbit interactions are known to have unusual 
effects on superconductivity 57,58 . The gate-tunable 
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where k B is the Boltzmann constant and T is the 
temperature. The total occupation number n = 
(1/1V 2 ) J2i,a,a ( C laa C iaa) is computed using the following 
relation 

n =±- 2 £ [\uLJ 2 f(En) + K a «\ 2 B-f(En))\ (7) 

n,i,a,cr 


FIG. 2. (Color online) The variation of the singlet and 
triplet pairing amplitudes (in meV), in intra-band and inter¬ 
band channels, with respect to the Rashba and atomic SOI 
strengths 7 (left column) and A so (right column) respectively. 
Other parameters are g = —0.6 eV, T = 0, g = 0.135 eV, 
h x 1 = 0.4 eV, h X 2 = 0.1 eV and rid = 0. 
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FIG. 3. (Color online) (a)-(d) The gap structures of the sin¬ 
glet and triplet superconductivity, in intra-band and inter¬ 
band channels. All pairing amplitudes are in meV. The max¬ 
imum (e) singlet and (f) triplet pairing amplitudes are plot¬ 
ted in the n — T plane. Other parameters are g = 0.135 eV, 
h x \ = 0.4 eV, h X 2 = 0.1 eV and rid = 0. 


FIG. 4. (Color online) (a)-(d) Variation of the pairing ampli¬ 
tudes with respect to external gate voltage V g . The maximum 
(e) singlet and (t) triplet pairing amplitudes are plotted in 
the V g — T plane. The white line in figure (e) is the experi¬ 
mental data, obtained from Ref. 28, showing the variation of 
the transition temperature T c with V g . Other parameters are 
g = 0.135 eV, h x i = 0.4 eV, h x 2 = 0.1 eV and rid = 0. 


Rashba SOI makes the interface q2DEG a potential can¬ 
didate for spintronic applications 59 as well as a play¬ 
ground for the search of non-trivial topological excitation 
such as the Majorana fermions 6 or the Skyrmions 60,61 or 
novel Hall phases 7 . We study the effects of Rashba SOI 
and atomic SOI on the intra-band and inter-band pairing 
as shown in FIG. 2. The pairing gap A a in the d xy or¬ 
bital is largely suppressed by the in-plane ferromagnetism 
while the pairing gaps A 5 and A c in the d yz , d zx orbitals 
dominate. As depicted in FIG. 2(a), (e), the pairing am¬ 
plitudes A b and A c decreases slowly with increasing both 
Rashba SOI strength 7 and atomic SOI strength A so be¬ 
cause the SOI enhances precession of electrons leading to 
slow reduction of the electron pairing. It is interesting 
to note that the superconductivity in the singlet inter¬ 
band channel and the triplet intra-band and inter-band 
channels is induced by the SOI when the pairing ampli¬ 
tudes At, or A c is finite as shown in FIG. 2(b)-(d),(f)-(h). 
This is because of the fact that the SOI breaks the spin- 
degeneracy and, in presence of the Zeeman splitting, the 
spin and orbitals nature of the electronic bands are re¬ 
organized and the pairings, in these channels, become 
energetically favourable. 

The temperature-variations of the pairing amplitudes 
are shown in FIG. 3(a)-(d). The pairing gaps reveal BCS 
nature in all the pairing channels except the triplet intra¬ 
band pairing A* which gets enhanced near the transition 
temperature. We plot the maximum pairing amplitudes 
in the singlet and triplet channels in the n—T phase-plane 
as described in FIG. 3(e)-(f). which show the supercon¬ 
ducting phases. It is evident that both the singlet and 
triplet pairing channels show a dome-shaped supercon¬ 


ducting phase and the triplet pairing becomes stronger 
towards higher values of career density. 

In the experiments, a gate voltage ( V g ) is tuned to con¬ 
trol the doping level of the q2DEG and a superconduc¬ 
tor to insulator transition with a quantum critical point 
V g = 140 V 28 is evinced by varying the gate voltage. It is, 
therefore, fascinating to study the variation of the pairing 
amplitudes in different channels with respect to V g and 
eventually to plot the phase-diagram in the V g — T plane. 
In appendix A, we derive the gate-voltage dependence of 
the career density and Rashba SOI. In FIG. 4(a)-(d), we 
plot the gate-voltage variation of the pairing amplitudes 
which reveal the superconducting transition at V g ~ 140 
V. It is interesting to note that the singlet intra-band 
pairing gaps and triplet intra-band and inter-band pair¬ 
ing gaps follow the nature of the variation of the Rashba 
SOI with respect to V g as shown in Fig. 8 (b). FIG. 4(e)- 
(f) depict the superconducting phases in the singlet and 
triplet channels respectively. The superconducting phase 
shown in FIG. 4(e) fits well with the experimental data 
except a quantitative mismatch of the transition temper¬ 
atures due to inherent overestimation problem of mean- 
field theory. The triplet pairing appears to begin at 
higher values in the gate voltage range. It has been found 
experimentally that the pairing gap A behaves differently 
than the transition temperature T c with respect to V g 50 . 
In the underdoped region, T c increases with V g while A 
decreases. This unusually different variation of A and 
T c with V g has been elucidated as due to the precise na¬ 
ture of the phonon spectral function a 2 g(ui ) in the un¬ 
derdoped regime and due to the entrance of new bands 
















































6 


(a) CT = 103 66 meV (b) = 103.66 mev 



FIG. 5. (Color online) The polar plots in figures (a)-(d) show 
the angular-variations of the maximum of singlet and triplet 
pairing amplitudes in the intra-band and inter-band channels 
respectively. The angle 9 is between the initial polarization 
direction due to the intrinsic ferromagnetism and the applied 
magnetic field. Figure (e) displays the variation of the super¬ 
conducting transition temperature T c with respect to the am¬ 
plitude B of the magnetic held. Parameters used are g, = —0.6 
eV, T = 0 eV, g = 0.135 eV, h x i = 0.4 eV, h X 2 = 0.1 eV and 
n d = 0. 


which suppress A in the overdoped regime 45 . 

Next we study the effect of external in-plane magnetic 
field on superconductivity. As shown in FIG. 5, the mag¬ 
netic field is applied in the interface plane along different 
directions (9) with respect to the initial direction of po¬ 
larization due to in-plane ferromagnetism and the pair¬ 
ing gaps are plotted in the polar plane of (A, 9). The 
singlet pairing amplitudes are increased along the oppo¬ 
site direction of that of the intrinsic polarization direc¬ 
tion while the triplet pairing amplitudes get strengthened 
along 9 ~ 60°, 300°. FIG. 5(e) plots the superconduct¬ 
ing transition temperature T c as a function of the ampli¬ 
tude B of the applied magnetic field. Remarkably, the 
enhancement of superconductivity by in-plane magnetic 
field has been reported experimentally in LAO/STO in¬ 
terface 51 . The enhancement of superconductivity is be¬ 
cause of the interplay between superconductivity, the in¬ 
trinsic ferromagnetism and the ferromagnetism induced 
by the applied magnetic field and implicates a possibility 



FIG. 6. (Color online) The variation of (a) the maximum 
pairing amplitude A max , (b) the gap-closing temperature T g 
and (c) the BCS ratio 2A rnax /T g with respect to the disorder 
strength W for different disorder concentrations n d . Other 
parameters are g = —0.6 eV, g = 0.135 eV, h x i = 0.4 eV and 
h X 2 = 0.1 eV. 


of a hidden superconducting phase above the transition 
temperature. The hidden superconducting phase is also 
inferred from the recent observation of the magnetic field 
assisted transient superconducting state observed in the 
interface q2DEG at 245 mK 52 . It is relevant to mention 
here that a hidden order has been proposed to be the 
precursor of superconductivity in Fe-pnictdes and high- 
Tc Cuprates 62 . It may, as well, be possible that a hid¬ 
den superconducting phase, arising from the competition 
between ferromagnetism and superconductivity, exist in 
LAO/STO hetero-interface. 


V. MONTE-CARLO STUDY OF THERMAL 
PHASE-FLUCTUATION 

The interface q2DEG is highly inhomogenous in nature 
and becomes superconducting at very low career concen¬ 
trations. Therefore, the superconductivity is prone to the 
detrimental effects of non-magnetic disorder and the role 
of phase-fluctuation becomes important. In a previous 
study 4 , we show that disorder, in fact, can help the antag¬ 
onistic ferromagnetism to live apart with superconductiv¬ 
ity in spatially phase-segregated regions at the interface. 
In FIG. 6, we show the variation of the maximum pairing 
gap A max , the gap-closing temperature T g and the BCS 
ratio 2A max /T g with respect to disorder strength W for 
different concentrations n d . With increasing W and n d , 
both A max and the BCS ratio decreases with fluctua¬ 
tion. On the other hand, T g fluctuates, with increasing 
W, within a narrow range of temperature and shows in¬ 
creasing tendency with increasing n d . It is evident that 
in the highly disorder limit, A and T g behave differently 
and, therefore, can no longer track the superconducting 
transition. It is required to go beyond the standard BCS 
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FIG. 7. (Color online) The variation of the phase correla¬ 
tion function Dij and the disordered and site-averaged pair¬ 
ing amplitude (|Ai|) with respect to (a) temperature T and 
(b) perpendicular magnetic field B z . Figure (c) depicts the 
density of states p(E) at different temperatures across the 
transition to the superconducting state, (d) The profile of 
the pairing amplitude (|Aj|) at T = 60 meV describing the 
localized Cooper pairs in the non-superconducting side of the 
transition (in figure (a)). A 14 x 14 lattice is used in the cal¬ 
culation. Other parameters are p = —0.6 eV, g = 0.135 eV, 
hxi = 0.4 eV, h X 2 = 0.1 eV, W = 0.8 eV and rid = 50. 


phase correlation function Dij = (cos (fa — fa)}, similar 
to that used in the XY-model to track the universal BKT 
transition, where i and j are sites separated by large dis¬ 
tance and ‘ ( ) ’ denotes the disorder and site averaged 
value. In FIG. 7(a)-(b), we show the variation of 
and (| Ai|) with respect to temperature T and perpendic¬ 
ular magnetic held B z respectively. In both the cases, 
the transition to the non-superconducting phase is dic¬ 
tated by a vanishing D ig and, on the other hand, (|Aj|) 
remains finite and fluctuates. The small finite value of 
(| A, |) in the non-superconducting phase is the signature 
of localized Cooper pairs 63 . The profile of the pairing gap 
in the two-dimensional space is shown in FIG. 7(d). The 
direct evidence of the presence of localized Cooper pairs 
in the insulating phase has been found in recent trans¬ 
port measurement 53 . FIG. 7(c) shows the single particle 
density of states, given by 


pW = n £ 


\<<xSH E - En) + K a afHE + E n ) 


n.i.Oi.cr 


at different temperatures across the superconducting 
transition. It is important to note that because of strong 
disorder, quasi-particle states, bound to the impurities, 
appear within the bulk superconducting gap and the sys¬ 
tem passes through a pseudo-gap like phase during the 
quantum phase transition 50 . 


prescription and incorporate the phase-fluctuation into 
account in the effective theory. 

We use a Monte Carlo (MC) technique in conjunction 
with self-consistent BdG formalism to study the ther¬ 
mal phase-fluctuations near the superconducting tran¬ 
sition driven by magnetic field or temperature. In the 
MC analysis, we consider only singlet intraband pairings 
in the d yz and d zx orbitals as they dominate over the 
pairings in other channels substantially and use a single 
pairing gap A = A£ ~ A s c to reduce the computation 
time. The complex superconducting order parameter is 
taken as A (n) = lAjle^, where |Aj| and fa are, respec¬ 
tively, the amplitude and phase of the pairing gap at site 
i. We start at high temperature with a specific disorder- 
configuration and reach low temperature upto 0.002 eV. 
At each temperature, the site-resolved pairing gaps |A*| 
are computed using the self-consistent BdG method, de¬ 
scribed in section III, and are fed into the MC update 
process which uses a free-energy minimization technique 
to determine the phases fa. In appendix B, we present 
the calculation of the free energy for a given BdG Hamil¬ 
tonian. 

In the weak to moderate disorder limit, the pairing gap 
parameter tracks the superconducting transition. How¬ 
ever, in the highly disordered limit, in which the interface 
q2DEG lies, the superconducting transition is indicated 
by the destruction of the global phase-coherence. To 
quantify the phase-coherence, we compute a long-ranged 


VI. CONCLUSIONS 

In this paper, we presented a three-band model for the 
superconductivity at the LaAlOa/SrTiOa (001) interface. 
We explored, using a multi-orbital BdG theory, the in¬ 
terplay between the superconductivity, ferromagnetism, 
and spin-orbit interactions. We also studied the role of 
thermal phase fluctuation using a Monte-Carlo method. 
The key findings of our analysis can be summarized as 
follows. 

The electron pairing in the d xy band is suppressed by 
the competing ferromagnetic order and the pairing in 
the d yz , d zx bands dominates. We find that the singlet 
pairing in the inter-band channel and the triplet pairing 
in both intra- and inter-band channel is induced by the 
SOI when the pairing amplitude in the singlet intra-band 
channel is finite. 

Gate-voltage has been found to affect supercocinduc- 
tivity at the interface quite strongly. We calculated the 
gate-voltage variation of the pairing amplitudes in differ¬ 
ent channels and extracted the singlet and triplet super¬ 
conducting phases in the V g — T plane and compared with 
experimental results with a fair degree of agreement. 

We observed an enhancement of superconductivity by 
external in-plane magnetic field. We found that the su¬ 
perconducting transition temperature depends on the di¬ 
rection of the applied in-plane magnetic field which is an 
expermentally verifiable prediction from our theory. The 
enhancement of superconductivity by applied magnetic 

















field is due to the interplay between ferromagnetism and 
superconductivity and suggests a ’hidden superconduct¬ 
ing phase’ above the transition temperature. The re¬ 
cent observation of magnetic field-assisted transient su¬ 
perconducting state at 245 mK 52 agrees with this hidden 
superconducting order and needs to be explored further 
experimentally 64 . 

Lastly, we studied the role of thermal phase fluctuation 
on the superconductivity using a Monte-Carlo method 
and found that there exist localized Cooper pairs in the 
non-superconducting phase beyond the quantum phase 
transition driven by perpendicular magnetic field. The 
same is also seen in the normal state just above the super¬ 
conducting transition temperature. The density of states 
reveal that, in the highly disordered situation, quasi¬ 
particle bound states appear within the bulk supercon¬ 
ducting gap and the system passes through a pseudo-gap 
like phase exactly as reported in the tunneling spectro¬ 
scopic measurement 50 . 
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Appendix A: Gate-voltage dependence 


The voltage V g , applied by back-gating to the STO 
substrate with the interface 2DEG grounded, controls 
the career density at the interface and the Rashba spin- 
orbit splitting 28,29 . The career density n changes with 
the applied electric field F according to the relation 

n{F) = — / e r {F)dF (Al) 

e J o 

where eo is the free-space permittivity, e r (F) is the field- 
dependent relative permittivity and e is the quantum of 
electronic charge. For bulk STO, e r (F) = 1/A(1 + ^F), 
where A = 4.097 x 10 ~ 5 and B = 4.907 x 10 “ 10 are 
temperature-dependent parameters, determined experi¬ 
mentally at 4.3 K by a first-order fit to the permit¬ 
tivity 65 . However, it has a singularity at F = —A/B 
which is circumvented, as prescribed in Ref. 66 , by tak¬ 
ing into account the second order term in the denom¬ 
inator. The corrected expression is given by e r (F) = 
l/A[Ci + C 2 fF + C 3 (f) 2 F 2 ], where C t (t = 1,2,3) are 
parameters to be found out by experimental fits. The 
career density can, therefore, be expressed as 


n{F) 



tan 1 


2C 3 bf+c 2 a 

A^fCi 



eBy/Ci 



FIG. 8. (Color online) Gate voltage modulation of (a) the 
career concentration and (b) the strength of Rashba spin- 
orbit splitting fitted with the experimental data 28,29 . The 
variation of (c) the total occupation number as a function of 
the chemical potential and (d) the chemical potential as a 
function of the gate voltage at zero temperature. A career 
concentration of 10 1 ' 5 cm -2 is equivalent to 0.016 el./u.c 67 . 


where C 4 = 4 C 1 C 3 — Cf. The gate-voltage V g enters in 
the above expression via the relation F = V g /d + Fq, 
where d = 0.5 mm is the thickness of the STO layer and 
F 0 = 1.2 x 10 5 V/m is the initial electric field at the 
interface in the absence of the gate-voltage. Comparison 
with the experimental data 28 yields C\ = 4.25, C 2 = 
—0.37 and C 3 = 0.29. Fig. 8 (a) describes the gate-voltage 
modulation of the career density. 

Furthermore, the gate voltage tunes the Rashba spin- 
orbit splitting at the interface 29 . In Ref. 67, the gate- 
voltage dependence of the Rashba spin-orbit interaction 
is derived using a Kane k • p approach. However, it is in¬ 
teresting to note that the DFT studies in Ref. 38 reveals 
that the Rashba spin-orbit interaction at the interface is 
very different from usual Rashba type spin-orbit inter¬ 
action found in semiconductor hetero-interfaces. Here, 
we fit the gate-voltage dependence of the Rashba split¬ 
ting amplitude observed in Ref. 29, to the equation 
7 = a + b(V g + 200) 2c within Fg-range [—200,200] V, 
with a = 0.0016, b = 3.6 x 10 -11 and c = 1.8 as plotted 
in Fig. 8 (b). 


Appendix B: Free energy of BdG Hamiltonian 


The free energy of an inhomogenous superconductor is 
given by 68 


(A2) 


*={X e/f )-TS 


(Bl) 
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where ( r K e ff) is the effective Hamiltonian of the super¬ 
conductor, T is the temperature and S is the entropy. 

The effective mean-field Hamiltonian of a single-orbital 
system is written as 

-H-eff = Hij c i<j c icr + y^[^ c jt c l|. + h.C ■] 4 jj— (B2) 
ij,u i 

where H,.j is the tight-binding Hamiltonian containing 
kinetic energy terms, Ai is the complex superconducting 
pairing gap at site i and U is the strength of the attractive 
electron-electron interaction. 

The effective Hamiltonian H e ff is diagonalized 
by using the Bogoliubov transformation Cj CT = 

En,*’ <4w W + V naa ,A /na' t0 ° btaln 

•^eff = Eg + yy ^n7io-'7n<T' (B3) 

n 

where E n are the eigen-values of the BdG equations and 


the ground state energy E g is given by 

E a = -Y, E n \<*\ 2 + E^f ( B4 ) 

i,n,cr i 

The entropy of an ideal gas of fermionic quasi-particles 
is expressed as 69 

§ = - £[/„ In fn + (1 - fn) ln(l - /„)] (B5) 

n 

Therefore, using (r na ,'y n <r') = fna’ and /„ T = f ni = 
f n , Eq. (Bl) becomes 

3=-Y J E n \v i na \ 2 + Y J ^+Y. E ^n (B6) 

i,n,(7 i n 

+ T y^[/n In f n + (1 — f n ) ln(l — /„)] (B7) 

n 

Evidently, the free energy can be calculated if the solu¬ 
tion of the BdG Hamiltonian is known. 
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